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THREE-DIMENSIONAL ELASTIC STRESS AND DISPLACEMENT ANALYSIS 
OF FINITE GEOMETRY SOLIDS CONTAINING CRACKS* 
by John P. Gyekenyesi and Alexander Mendelson 
Lewis Research Center 

SUMMARY 

The line method of analysis is applied to the Navier-Cauchy-equations of 
elastic equilibrium to calculate the displacement distributions in various 
bodies containing cracks. The application of this method to these equations 
leads to coupled sets of simultaneous ordinary differential equations whose 
solutions are obtained along sets of lines in a discretized region. When de- 
coupling the equations and their boundary conditioiis is not possible, the use . 
of a successive approximation procedure permits the analytical solution of . 
the resulting ordinary differential equations. The results obtained show 
considerable potential for using this method in the three-dimenaibnal ana- 
lysis of finite geometry solids and suggest a possible, extension of this tech- 
nique to nonlinear material behavior, 

INTRODUCTION 

Considerable progress has been made in recent years in the stress 
analysis of bodies containing flaws or cracks. However, most of the work 
on this subject has been based on the plane theory ©f elasticity (ref, 1), 

For a fracture specimen of finite geonietry the stress and displacement fields 
are highly three -diniensional,,' Hence, solutions of the general elasticity 

* Submitted to the International Journal of Fracture Mechanics., 



2 


field equations must be obtained when more reliable results are needed. 

Because of the geometric singularity associated with any crack type 
problem, there is almost no possibility of a simple closed form type of so- 
lution, For this reason, three-dimensional elastic solutions have been 
obtained only for a restricted class of problems (refs, 2 to 4), Recently, 
with the availability of large digital computers, a variety of numerical 
methods appeared in the literature; however, most of these methods have 
yielded only partial results. Among these approxiniate methods are the 
finite difference (ref, 5), the finite element (ref, 6), the 4irect potential ( 
(re"f, .7), the eigenfunction expansion (ref, 8), and the line method of ana- 
lysis (refs. 9: to ll). Of all these sblutiont technique?, wthejline method^of 
analysis appears to yield the most complete and accurate results for three- 

I 

dimensional elasticity problems. 

Although the concept of the line method for solving partial differential 
equations is not new (ref, 12), its application in the past has been limited to 
simple examples (ref, 13). The line method lies midway between completely 
analytical and discrete methods. The basis of this technique is the substi- 
tution of finite differences for the derivatives with respect to all the inde- 
pendent variables except one for which the derivatives are retained. This 
approach replaces a given partial differential equation with a system of si- 
multaneous ordinary differential equations whose solutions can then be ob- 
tained in closed form. These equations describe the dependent variable 
along lines which are parallel to the coordinate in whose direction the deriv- 
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atives were retained. Application of the line method is most useful when 
the resulting ordinary differential equations are linear and have constant 
coefficients. 

An inherent advantage of the line method over other numerical methods 
is that good results are obtained from the use of relatively coarse grids. 

This use of a coarse grid is permissible because parts of the solutions are 
obtained in terms of continuous functions. Additional accuracy in normal 
stress distributions is derived from the fact that they are expressed as first- 
order derivatives of the displacements and these derivatives can be analyti- 
cally evaluated. Inherently inaccurate numerical differentiation is required 
only for evaluating the shear stresses, but this presents ho important loss of 
accuracy since they are an order of magnitude smaller than the normal 
stresses. For problems with geometric singularities, additional accuracy 
is derived from using a displacement formulation since the resulting defor- 
mations are not singular. 

It is the purpose of this report to present a simple and systematic 
approach to the elastic analysis of three-dimensional, finite geometry solids 
containing traction-free cracks. The need for these specific solutions has 
existed for a number of years in fracture toughness testing. Problems that 
are most conveniently described in rectangular Cartesian coordinates as 
well as circular geometry solids are treated. 
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:^EDUCTIONOF THE NAVIER-CAUCHY EQUATIONS TO SYSTEMS 
■ OF ORDINARY DIFFERENTIAL EQUATIONS 
Within the framework of linearized elasticity theory, the equations of 
elastic equilibrium in terms of displacements are 


(X + G) — + GV^U = 0 

( 1 ) 

9x 


(X +G) — + GV^,v = 0 

( 2 ) 



(X + G)- de ^ ^ 0 

( 3 ) 

3z 



where the body forces are assumed to be zero^and the dilatation is 


9u 9v 9w 
€=—+ — + — 


( 4 ) 


9x 9y 9z 

For a finite geometry solid with rectangular boundaries, we construct 
three sets of parallel lines (fig, 1(a)). Each set of lines is parallel to one 
of the coordinate axes and thus perpendicular to the corresponding coor- 
dinate plane. An approximate solution of equation (1) can then be obtained by 
developing solutions of ordinary differential equations along the x-directional 
lines. As seen in the figure, there are a total of f = NY x NZ such lines 
where NY is the number of lines along the y-direction and NZ is the number 
of lines along the z-dire;ction in a given plane, respectively. We define the 
displacements along these lines as U 2 , . , , , u^. The derivatives of the 
y-difectional displacements on these lines with respect to y are defined as 


HI 


> V 


, v’ ^ , and the derivatives of the z-directional displacements 
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with respect to z are defined as w*|j, ^’|2’ ■ ’ * > These displace- 

ments and derivatives can then be regarded as functions of x only since 
they are variable^ on x-directional lines. When these definitions are used, 
the ordinary differential equation along a generic line ij (a double subscript 
is used here for simplicity of writing) in figure 1(b) may be written as 


^2 

“ij , a - 2v) 


2 2 \ 1 / ■ ■ V 

' y “ z ' 


“z 


f..(x) 

-ii— =0 


2(1 - u) 


( 5 ) 


where 


fij(x) = 


dv* 

dx 


dw' 


1] 


dx 


ij 


( 6 ) 


v’ 


d*ir 


and 


dy 


dz 

Similar differential equations are obtained along the other x-directional 
lines. Since each equation has the terms of the displacements on the sur:-? 
rounding lines, these equations constitute a system of ordinary differential 
equations for the displacements Uj^, Ug, . . . , u^. 

The set of S. second order differential equations represented by (5) 
can be reduced to a set of 2it first order differential equations by treatii^ 
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the (derivatives of the u’s as an additional set of £ unknowns, i,e, defining 

dUi dUg 

U/, I, = — f U- »= - 3 — , etc. 

£+1 dx ’ £+2 dx' 

The resulting 2£ equations can now be written as a single first order 
matrix differential equation 


( 7 ) 


^ = AjU+-Rj(x) (8) 

where U and R are column matrices of 2£ elements each and 

is a 2£x2£ matrix of the constant coefficients appearing in equations (5) 

and (7). 

In a similar manner, to solve equations (2) and (3) ordinary differen- 
tial equations are constructed along the y- and z-directional lines respec- 
tively, TheTse equations are also expressed in an analogous form to equa- 
tions (8); they are 

^ = A2V,S(y) ^ (9) 

^=A3W + T(z) (10) 

Equations (8) to (10) are linear first-order ordinary matrix differential 
equations. They are, however, not independent, but are coupled through 
the vectors R, S and T whose components are given by equations similar 
to (6). The elements of the coefficient matrices Aj^, Ag and A^ are all 
constants, being functions of the coordinate increments and Poisson's ratio 
only. 
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]Sroting that a second-order differential equation can satisfy only a total 
of two boundary conditions and since three-dimensional elasticity problems 
have three boundary conditions at every point of the bounding surface, some 
of the boundary ddta must be incorporated into the surface line differential 
equations. Hence, conditions of normal stress and displacement are en- 
forced through the constants of the homogeneous solutions while shear stress 
bounflary data must be incorporated into the differential equations of the 
surface lines. The application of the specified shear conditions permits the 
use of central difference approximations when surface line differential eq- 
uations are constructed. The details of constructing these equations are 
found in reference 

EQUATIONS IN CYLINDRICAL COORDINATES 

! 

The Navier- Cauchy equations in cylindrical coordinates are written as 



where the body forces are assumed to be zero and the dilatatiomis 

= ^ + + ^ ( 12 ) 

9r r 36 r 3z 

For a finite geometry body with circular boundaries, we construct 
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three sets of parallel lines in the direction of the coordinates as shown in 
figure 2. Approximate solutions of equations (11) can then be obtained by 
developing solutions of ordinary differential equations along the radial, 
circumferential, and axial lines, respectively. For the first of equations 
(11), we define the displacements along the radial lines as u^, . . • 

The derivatives of the circumferential displacements on these lines with re- 
spect to 0 are defined as v’ v' , v'j^ , and the derivatives of the 

axial displacements with respect to z are defined as w’ w* • • • > 
w'jg . These displacements and derivatives can then be regarded as functions 
of the radius only, since they are variables on radial lines. If these def- 
initions are used, the ordinary differential equation along a generic, radial 
line ij of figure 2 may be written as 


d^u.. ^ du.. u.. 

1 ] 1 1 , 1] 1.1 

, (1 - 2l>) 

,2 r dr 2 

dr r 

2(1 - P) 


2 2 \ 1 , 

+ — u.. + (u. ..^ . :+ u. .) 

2 i ,2 .,2 / 2 . 2 

*>8 V '■ “’0 


+ ((u. . . + u. . J 

2 k]+l 


f..(r) 

1] ' 

2(1 - v) 


= 0 (13) 


where 



(4i/ - 3) , 

1 dv' 

1 . 

dw’ 

1 

V 

2 

r 

♦ « 1 

r dr 



dv 


v’ = 

: 



de 



( 14 ) 


and 



9 


dz 

Similar differential equations are obtained along the other radial lines. 

Since each equation has the terms of the displacements on the surrounding 
lines, these equations constitute a system of ordinary differential equations 
for the displacements u^, Ug, - . . , u^. 

Noting that a second-order differential equation can satisfy only a total 
of two boundary conditions, the shear stress boundary data must again be 
incorporated into the surface line differential equations. For the first radial 
line of figure 2, the use of zero shear stress boundary conditions in the 
radial direction on the r, z and r,9 coordinate planes gives, respectively, 
the following imaginary radial line displacements: 


ou dVi 

"lee = “2 2he>- ^ 




- VI 1 


=u„o,i+2h 

10n N0+1 z , ' 

dr 


(15) 


Equations (15) must then be used in the application of central difference 
approximations when the ordinary differential equation for the first radial 
line is generated. Additional details on the construction of these equations . 
can be found in reference 9. 

As before the system of jf second order equations given by equation (13) 
can be written as a single first order matrix equation 


^ = A (r)U +R(r) 
dr ^ 


(16) 
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where in this case the elements of U are defined by 


Uj = ruj, U2 = rUg, = ru^ 


d(ruJ 


U 


f+1 




r dr 


U 


f+2 


I d(rug)^ 

r dr 


• • • 


U 


2jf 


^ d(rUjj) 

r dr 


(17) 


and the elements of the coefficient matrix A^(r) are no longer constant, 
but are functions of r. The corresponding equations for V and W are 
the same as (9) and (10) with y replaced by 0, the coefficient matrices 
being constant. 

SOLUTION OF THF SYSTEMS OF ORDINARY DIFFERENTIAL EQUATIONS 


The. systems of ordinary differential equations (8) to (10) and (16) can be 

solved by any of a number of standard techniques. The method used herein 

was basically the matrizant or Peano-Baker method of integration (ref. 14). 

For equation (8) the solution can be written as 

A-x A.x /^'x -A 77 

U(x) = e U(0) + ® ® (1^) 

with similar solutions for equations (9) and (lO). U(0) is the initial value 
vector, determined from the boundary conditions. The conversion of given 
boundary data into required initial values is discussed in more detail in 
reference 9. 

Ax 

The matrizant, el is generally evaluated by its matrix series. For 
larger values of x, when convergence becon^qs slow, additive formulas 
may be used. Iii addition similarity transformations can be used to diago- 
nalize the matrix Aj. ' These various techniques for improving the accuracy 
are discussed in detail in reference 9. 
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For the case of equation (16) the solution is more complicated, since 
we no longer have constant coefficients. The solution may be expressed 
as follows (ref. 14); 

U(r) = J2(A^)U(r^) + P S2(A^)j R:(T?)d77 (19) 

J^o 

where r^ is the initial value of r which may or may not be zero. The 
matrizant J^(A^) is given by the infinite matrix integral series 

fi(Ap = I+ r A^ir)^)dr]^+ P A'^(n2)^V2 P 

Jo Jo Jo 


A^(T73)dr?3 










which of course becomes very difficult to evaluate. However, by substi- 
tuting the matrix A^ into equation (20), it can be seen by inspection, (ref. 
9), that n(A^) can be partitioned into four submatrices which satisfy the 
following equations; 
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and K is the sub matrix of A obtained from the identity:; 
i r r 



Instead of evaluating the matrizant n(A^) by means of equation (20), 
it can be evaluated by solving equations (21), This was done for the examples 


presented herein, using the single step Runge-Kutta integration method for 
solving equations (21). 

Since equation;^ (8) to (10) and their boundary conditions are highly 
coupled, it is generally impossible to directly evaluate their solutions. 
Thus, a successive approximation procedure must be employed where as- 
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sumed values must be used initially for the required unknowns. The cyclic 
re substitution of the obtained solutions into the coupling vectors and the 
boundary conditions will usually converge to the correct solution, depending 
mainly on the accuracy to which the required matrizant can be evaluated. 

Once the displacement field in the body has been calculated and the 
successive approximation procedure has converged, the normal stress 
distributions can be obtained directly by using the stress -displacement 
equations. The shear stresses, however, can be evaluated only through 
finite difference approximations for the required displacement gradients. 

NUMERICAL RESULTS 

Solid Cylindrical Bar with Penny-Shaped Crack 
Figure 3 shows a cylindrical bar containing a penny-shaped crack and 
loaded by a uniform normal stress distribution. For problems with axisym- 
metric geometry, the circumferential displacement is inherently zero at 
every point and all the remaining variables are independent of the circum- 
ferential coordinate 0, The two sets of parallel lines needed for the soluir 
tion of this problem are also shown in this figure. Note that the crack edge 
is assumed to be midway between adjacent nodes, specifying normal stress 
and displacement boundary conditions, respectively. 

The solution of this problem was obtained by using two different sets 
of lines along the coordinate axes so that the convergence of the finite differ- 
ence approximations could be checked. Selected results are shown in figures 
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4-6. Figure 4 shows that a relatively course grid of nine axial lines and 
nine radial lines gave almost identical results to those obtained using a 
16 by 13 grid. Figures 5 and 6 compare the results obtained herein with 
analytical results of references 15 and 16 for infinite bodies and infinitely 
long finite radius cylinders with cracks. Both figures show that as the 
dimensions of the cylinder increase relative to the crack radius, the results, 
from the method of lines approach the analytical results, as expected. 

Table I shows the axial stress distribution for one case. The accuracy 
of the normal stress boundary condition)^ can be noted from the listed results. 

Annular Plate with Internal Surface Cracks 

A much more difficult problem of a thick annular plate with four sym- 
metrically situated internal surface cracks and loaded by a uniform radial 
tensile load, is shown in figures 7 and 8, Because of symmetry only one- 
sixteenth of the orginal plate has to be discretized. Only a relatively 
coarse grid of 4 by 4 by 4 lines was used. Some of the results are shown 
in Tables 2 to 5, 

Because of the use of a coarse grid, the numerical results for this 
example are somewhat inaccurate in magnitude but they do indicate some 
previously unknown variations in the stress field for this problem. This 
conclusion is possible in that the line method does not usually require a fine 
grid for good results as was shown in the previous section. These results 
also demonstrate that the method of lines permits the computation of the 
displacement and stress fields for a general three-dimensional problem. 
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APPjLICATION TO TENSILE FRACTURE SPECIMENS CONTAINING CRACKS 
A great amount of experimental work has been done in fracture mechan- 
ics (ref. 17) through the use of crack-notched specimens. In the past, 
many different types of specimens have been used to determine a material’s 
fracture toughness. The most common early specimens employed in these 
tests were the center-cracked and double -edge -notched bar specimens. 
Figures 9(a) and 10(a) show the finite rectangular bars with through -thick- 
ness, traction-free central and double-edge cracks, respectively. Be- 
cause of the symmetric geometry and loading, only one-eighth of the bars 
has to be discretized as shown in figures 9(b) and 10(b). 

NUMERICAL RESULTS 
Center- Cracked Tensile Fracture Specimen 
The solution of this problem was obtained by using two different sets 
of lines along the coordinate axes so that the convergence of the finite dif> 
ference approximations could be checked, In a given direction, uniform 
line spacing was used in all computations with no other restriction being 
placed on the selection of the grid size. The crack edge location with re- 
spect to the imposed grid was assumed to be halfway between nodes spec- 
ifying normal stress a^nd displacement boundary conditions, respectively. 

The successive approximation procedure required for decoupling the three 
sets of ordinary differential equations was terminated when the difference 

between successively calculated displacements at every point was less than 

“6 

a preset value (10 ). 
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Selected results of the dimensionless displacements are listed in tables 6 
to 8. The dimensionless crack opening displacement;isrshown.iri fi^re 11. - 
Inspection of figures 11(a) and (b) shows that the grids of lines have been 
sufficiently refined when the results of the 64 by 112 by 112 grid were calcu- 
lated. Figure 11(a) also contains the results of the plane elasticity solutions 
obtained by Mendelson, Gross, and Srawley (ref. 18). It is noteworthy that 
the results correspond to elliptical crack profiles in all cases. As can be 
noted from figure 11(c), the finite length of the bar has a very noticeable 
effect on the crack opening displacements for values of L<3a. Figure 11(d) 
shows the variation of the crack opening displacements across the thickness 
of the bar. 

The dimensionless normal stress distributions in the crack plane are 
shown in figures 12 to 14 as a function of 'both: thie x- and z-coor-dirtates. • 

The results in these figures clearly indicate the singular nature of the nor- 
mal stresses near the crack edge. As expected, the stress normal to the 
crack plane increases most rapidly near the crack front. A plot of these 
stresses in the z -direction also indicates a central region of uniform stress 
and a boundary layer through which the stresses decrease to the surface 
values. These same results show that as x increases the stress field ap- 
proaches a uniaxial state of stress which indicates that the cause of a tri- 
axial stress field in the bar is the through -thickness central crack. 
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Double-Edge Crack Tensile Fracture Specimen 
The crack opening displacements for the problem of figure 10(a) are 
presented in figure 15<. Comparing figures 15(a) and (b) shows that the finite 
difference approximations have sufficiently converged when the results of 
the finer grid were calculated. Contrary to the central crack problem, the 
crack opening displacements in figure 15 are independent of the z-coordinate. 
Similar results were obtained by Cruse and Van Buren (ref, 7) for the single- 
edge-crack bar specimen. 

STRESS INTENSITY FACTOR 

It is customary in fracture mechanics to describe the plane elasticity 
crack opening displacement as a superposition of three basic deformation 
modes (ref. 1). Since the problems shown in figures 9(a) and 10(a) have geo- 
metric symmetry and are symmetrically loaded^ , only the opening mode of 
crack displacement is obtained. In terms of the stress intensity factor for 
the opening mode Kj, the plane elasticity crack displacements near the 
crack tip are given by (ref. 1) 

plane strain (22) 

plane stress (23) 

where R here is the distance from the crack edge. Since three-dimension- 
al problems are neither in a state of plane strain nor in a atate of plane 
stress, the definition of a stress intensity factor for these problems must 
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be first established. Note that by definition the plane stress and plane 
strain stress intensity factors are equal while the displacements are approx- 
imately 12.5 percent different for v = 1/3. Since the results indicate that 
for a sufficiently long bar most of the bar: in‘ the thickness direction is ap- 
proximately in a state of plane strain, equation (22) is selected to calculate 
the stress intensity factor. Rearranging this equation so that the dimen- 
sionless crack opening displacements can be used leads to 


where 



(24) 


4(1 - 

r~ 

Ev 27ra 


A plot of equation (24) as vH/a — 0 can then be used to calculate Kj, 

Since the crack opening displacement is a function of the thickness variable, 
the previously defined stress intensity factor varies in the z-direction. 

It should be noted that this description of Kj is completely arbitrary and 
that it is questionable if it has any real signMicance in three-dimensional 
elasticity problems. However, the values of Kj. are presented here so 
that a comparison is possible with published plane strain solutions (ref, 1). 
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Figure 16 shows the variation across the thickness of the stress in- 
tensity factor for the center -era eked bar. Note that Kj is a maximum, 
at the surface and a minimum near the center. The value of Kj for the 
double -edge -crack specimen is independent of the z coordinate since the 
crack opening displacement is constant across the thickness. 

Although the stress intensity factors for these problems could be de- 
termined with reasonable accuracy, the associated type of singularities 
are difficult to evaluate because values of the normal stresses are needed 
within a distance of 0. 05a or less from the crack edge. With the equal 
spacing of lines used in these examples, the minimum node location for these 
problems is about 0. 06a, For this range of crack edge distance R, the 
singularity of the stresses is not defined. 

CONCLUSIONS 

The line method of analysis presented affords a practical straight for- 
ward way for analysis of three-dimensional crack problems, at least for 
bodies with reasonably regular boundaries. Because parts of the solution 
are obtained as continuous functions along the lines chosen, relatively good 
accuracy can be obtained with coarse grids. In addition it should be noted ' 
that the introduction of plasticity into the analysis could be accomplished 
by merely changing the coupling terms appearing in equations (14) to (16). 
Since these have to be determined by an iterative process in any case, it 
would seem possible to solve the elastoplastic problem by a simple exten- 
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Sion of the present method. Whether this approach is practical requires 
further investigation. 
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TABLE 1. - NONDIMENSIONAUZED AXIAL STRESS FOR 

SOLID CYLINDRICAL BAR WITH PENNY-SHAPED CRACK 
UNDER UNIFORM NORMAL TENSION 


|_a = 1. 0, b = 1. 77, L =1.68 (16 axial and 13 radial lines)] 


z 

r 

■ 

■ 


0.000 

0.235 



0.941 

1.060 

— 

1.294 



0.00 

0.000 

0.000 

O.QOO 

0.000 

0.000 


1.512 



.28 

.095 

.093 

. 146 

.319 

.874 

1.513 

1.365 


[H[l^ 

.56 

.277 

.295 

.386 

.592 

.950 


1.230 

1. 186 

1.172 

.84 

.515 

.542 

.624 

,770 

.957 

1.041 

1.122 

1.136 

1. 159 

1. 12 

.736 

.759 

.808 

.885 

.972 

1.010 

1.059 

1.082 

1.092 


.900 

.915 

.933 

.960 


1.003 

1.022 

1.036 

1.037 


.999 

.998 

.999 

.999 

1.000 

.999 

.995 

.993 

.990 


















TABLE 2. - DIMENSIONLESS ORCUMFERENTIAL 
DISPLACEMENTS Ev/a^b FOR ANNULAR PLATE 
WITH INTERNAL SURFACE CFiACKS UNDER 
UNIFORM RADIAL TENSION ON 


OUTSIDE SURFACE 



















TABLE 3. - DIMENSIONLESS RADIAL STRESS 
DISTRIBUTION FOR ANNULAR 

PLATE WITH INTERNAL SURFACE 
CRACKS UNDER UNIFORM 
RADUL TENSION ON 


OUTSIDE SURFACE 


























TABLE 4. - DIMENSIONLESS aRCUMFERENTIAL 
STRESS FOR ANNULAR PLATE WITH 

INTERNAL SURFACE CRACKS UNDER 
UNIFORM RADIAL TENSION ON 


OUTSIDE SURFACE 


















TABLE 5. - DIMENSIONLESS AXIAL STRESS 
FOR ANNULAR PLATE WITH INTERNAL SURFACE 
CRACKS UNDER UNIFORM RADIAL TENSION 
ON OUTSIDE SURFACE 





TABLE 6. - DIMENSIONLESS x-DIRECTIONAL DISPLACE- 


MENTS Eu/a a FOR RECTANGULAR BAR UNDER 
o 

UNIFORM TENSION CONTAINING THROUGH- 
THICKNESS CENTRAL CRACK 

[a = 1. 0; b = 2. 0; L = 1. 75; t = 1. 5 (64 by 112 by 112 x-, y-, 
z- directional lines, respectively)^] 


y 

X 

z 

0. 

00 

0.31 

0. 77 

1.23 

1.69 

2.00 

x-Directional displacements, Eu/or^a 

0. 00 

0. 

00 

-0.478 

-1.047 

-1.320 

-1.391 

-1.437 

0. 

00 

.50 



-.123 

-.307 

-.678 

-.999 

-1.103 



1.00 



-.002 

-.059 

-.245 

-.463 

-.572 



1.75 



.220 

.447 

.489 

.404 

.324 



o.ooi 

0. 

00 

-0.494 

-1. 084 

-1.305 

-1.331 

-1.383 

1. 

07 

.50 



-.111 

-.282 

-.636 

-.949 

-1.061 



1.00 



.011 

-.034 

-.213 

-.431 

-.546 



1.75 



.241 

.489 

.534 

.438 

.350 



0.00 

0. 

00 

-0.410 ' 

n 

-1.196 

1 

-1.324 

-1.340 

1. 

50 

.50 

1 

1 


-.043 


-.556 

-.934 

-1.061 

j 


1.00 

1 


.036 

.013 

-.183 

-.424 

-.546 



1.75 

1 

1 


.238 

.483 

.525 

.428 

^ 1 

.338 

j 

































TABLE 7. - DIMENSIONLESS y-DIRECTIONAL DISPLACEMENTS 
Ev/a^a FOR RECTANGULAR BAR UNDER UNIFORM 
TENSION CONTAINING THROUGH>THICKNESS 
CENTRAL CRACK 


(a = 1.0; b = 2.0; L = 1.75; t = 1.5 (64 by 112 by 112 x-, y-, 
z-directional Lines, respectively )[] 


z 

y 


X 


0.00 

0.25 

0.50 

1.00 

1.50 

1.75 



y-Directional displacements, Ev/ct^ 

a 

0, 00 
.43 
1,07 
1.50 

2.742 

2.769 

2.946 

3.258 

2.892 
2. 912 
3.070 
3.343 

2.997 

3.018 

3.157 

3.397 

3.227 

3.246 

3.369 

3.569 

3.536 
3. 556 
3.678 
3.871 

3.688 

3.708 

3.835 

4.034 

0. 00 

1 

0.00 

.43 

1.07 

1.50 

1.703 

1.718 

1.812 

1,974 

1.849 

1.863 

1.952 

2.088 

2.055 

2.069 

2.158 

2.283 

1 

2.526 

2.542 

2.631 

2.752 

2.972 

2.988 

3.080 

3.205 

3.173 

3.189 

3.282 

3.410 

0.77 

1 

0. 00 
.43 
1,07 
1.50 

0. 000 

) 

0.328 

.332 

.348 

.343 

0. 689 
.695 
.722 
.722 1 

1.362 

1.371 

1.409 

1.432 

1.924 
1.932 
1.973 
2. 007 

2.178 

2.186 

2.225 

2.257 

1.54 

1 



Bi 

0. 328 
.333 
.353 
.353 

0.832 

.829 

.860 

.874 

1.328 

1.334 

1.362 

1,383 


1 







TABLE 8. - DIMENSIONLESS z-DIRECTIONAL DISPLACEMENTS 

Ew/a a FOR RECTANGULAR BAR UNDER UNIFORM 
0 

TENSION CONTAINING THROUGH -THICKNESS 
CENTRAL CRACK 


[a = 1. 0; b = 2. 0; L = 1. 75; t i 1. 5 (64 by 112 by 112 x-, y-, 
z-directional lines, respectively )] 


X 

z 

y 



0.00 

0.43 

0.64 

1.07 

1.28 

1.50 





z-Directional displacements, Ew/a^a 



0.00 

o.ool 

0.036 

0. 068 

0.185 

0. 272 

0.366 

0. 

0 

.77 



-. 055 

-.081 

-.123 

-.124 

-.095 



1.54 



-.161 

-.249 

-.459 

-.581 

-.704 



2.00 

1 


-.179 

-.270 

-.441 

-.512 

-.568 


0.00 

0.00 

-0, 042 

-0. 059 

-0. 077 

-0. 082 

-0. 090 

0. 

5 

.77 



-.097 

-.147 

-.247 

-.296 

-.342 



1.54 



-.162 

-.247 

-.429 

-.522 

-.610 



2.00 



-.176 

-.264 

-.430 

-.505 

-.574 

• t 

0.00 

0. 00 

-0. 123 

-0.184 

-0.307 

-0.370 

-0.435 

1. 

0 

.77 



-.144 

-.217 

-.366 

-.442 

-. 517 



1.54 



-.168 

-.253 

-.423 

-.507 

-.590 



2.00 

1 


-.176 

-.262 

-.425 

-.501 

-.577 



0.00 

0.00 

-0.222 

-0.339 

-0. 595 

-0. 738 

-0.886 

1. 

75 

.77 



-.212 

-.320 

-.546 

-.662 

-.780 



1.54 



-.191 

-.285 

-.465 

-.549 

-.629 



2.00 



-.190 
.1 

-.281 

-.446 

-.521 

-.593 







(a) Three sets oMines parallel to X- y- and z -coord mates and 
perpendicular to corresponding coordinate planes. 



(b^ Set of interior lines parallel to x-coordinate. 


Figure 1. - Sets of lines parallel to Cartesian coordinates. 
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Figure 2 . - Sets of lines is direction of cylindrical coordinates. 




Figure 3. “ Solid cylindrical bar with penny-shaped crack. 
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ib) Dimensionless axial displacement distribution (9 by 9 grid). 

Figure 4. - Dimensionless axial displacement distribution 
for solid cylindrical bar with penny-shaped crack. 
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Figure 5, - Dimensionless crack opening displacements for 
solid cylindrical bars with penny-shaped cracks of various 
lengths and radii. 






Figure 6. - Dimensionless maximum crack opening as 
function of crack to cylinder radius ratio. 



Figure 7. - Annular plate with interna! surface cracks under 
uniform external tension. 
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la) Rectangular bar wfth throuah ‘thickness double edge cracks. 



(b) Discretized region of rectangular bar with double-edge cracks. 

Figure 10. - Rectangular bar with through -thickness double-edge 
cracks under uniform tension. 



(b) Dimensionless crack opening displacement 
05x70x98 grid). 

Figure 11, - Dimensionless crack opening displacement 
for a rectangular bar under uniform tension con- 
taining a through -thickness central crack. 



a = 1.0 
b = 2,0 



0 LO 2.0 3.0 4.0 

L 


(cl Variation of crack opening displacement 
with bar length. 
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(d) Dimensionless crack opening displace* 
ment variation across bar thickness. 


Figure 11. - Concluded. 



(a) Dimensionless X -directional normal stress 
variation along bar width. 
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(bl Dimensionless x-dlrectional normal stress 
variation across bar thickness. 

Figure 12. - Dimensionless x -directional normal 
stress distribution in the crack plane for a 
rectangular bar under uniform tension con- 
taining a through -thickness central crack. 
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(a) Dimensionless y -directional normal stress 
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(b) Dimensionless y-directional normal stress 
variation across bar thickness. 

Figure 13. - Dimensionless y-directional normal 
stress distribution in the crack plane for a 
rectangular bar under uniform tension con- 
taining a through -thickness central crack. 



o|y*0 



(a) Dimensionless z-directional normal stress 
variation across bar width. 



(b) Dimensionless z-dlrectional normal stress 
variation across bar thickness. 


Figure 14. - Dimensionless z-directional normal 
stress distribution in the crack plane for a 
rectangular bar under uniform tension con- 
taining a through -thickness central crack. 
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h. la) Dimensionless crack opening displacement 
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(b) Dimensionless crack opening displacement 
(35x70x98 grid). 

Figure 15. - Dimensionless crack opening 
displacement for a rectangular bar under 
u nifor m te ns ion conta i n i ng through ' 
thickness double edge cracks. 
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Figure 16. - Variation of stress intensity factor 
Kj across thickness for rectangular bar 
under uniform tension containing through- 
thickness central crack. 



